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Abstract 

Time-dependent density-functional theory (TDDFT) is widely used to describe electronic excita- 
tions in complex finite systems with large numbers of atoms, such as biomolecules and nanocrystals. 
The first part of this paper will give a simple and pedagogical explanation, using a two-level system, 
which shows how the basic TDDFT formalism for excitation energies works. There is currently 
an intense effort underway to develop TDDFT methodologies for the charge and spin dynamics 
in extended systems, to calculate optical properties of bulk and nanostructured materials, and 
to study transport through molecular junctions. The second part of this paper highlights some 
challenges and recent advances of TDDFT in these areas. Two examples are discussed: excitonic 
effects in insulators and intersubband plasmon excitations in doped semiconductor quantum wells. 

Keywords: time-dependent density-functional theory; excitation energies; semiconductors; excitons; quan- 
tum wells; intersubband plasmons 
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I. INTRODUCTION 



Many important areas of experimental and theoretical physics, chemistry, and materials 
science require an understanding of the electronic excitations of atomic or molecular systems, 
nanostructures, mesoscopic systems, or bulk materials.-^ 1 ^ A wide variety of spectroscopic 
techniques is being used to characterize the electronic structure and dynamics of these 
systems by probing their excitation spectra. The performance of any nanoelectronic device, 
such as a molecular junction, is dominated by its electronic excitations.- 1 ^ 

Time-dependent density-functional theory (TDDFT) 7 is an increasingly popular, uni- 
versal approach to electronic dynamics and excitations. Just like ground-state density- 
functional theory (DFT), which is based on a set of rigorous theorems^ proving a one-to-one 
correspondence between ground-state densities and potentials, there is a similar existence 
theorem for TDDFT, due to Runge and Gross,— which establishes the time-dependent den- 
sity as a fundamental variable. 

The usage of TDDFT as a practical tool to calculate excitation energies started in the 
mid-90's with the groundbreaking work of Petersilka et a/.— and Casida.— In the Casida- 
formalism, the excitation energies are obtained from the following eigenvalue problem: 

(1) 

) 1 / \ Y I 

where the matrices A and K are defined as follows: 
and 

Kiaa,i'a'a' = J d z r J d 3 r'ip* a (r)(f aa (r) 





/I /xc, 



r — r 



Vi'v' A r ')- ( 3 ) 



Here, ip,- ia and are the Kohn-Sham orbitals and eigenvalues coming from a self-consistent 
ground-state DFT calculation; we use the standard convention that i, i', . . . are indices of 
occupied orbitals and a, a', . . . refer to unoccupied orbitals. / XCj(Jcr /(r, r', cu) is the so-called 
exchange-correlation (XC) kernel,— which is in general a frequency-dependent quantity, but 
in practice is often approximated using frequency-independent expressions. This is known 
as the adiabatic approximation. 
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The eigenvalues to of Equation ([I]) are, in principle, the exact excitation energies of the 
system, provided one would start from an exact Kohn-Sham ground-state calculation and 
then use the exact / xc ,o-cr'- In practice, of course, static and dynamical XC functionals need 
to be approximated. The formalism of Equations ([I])-© can also be recast in the shape of 
an eigenvalue problem for the squares of the excitation energies: 



where u iaa = e a(7 - e ia . 

This TDDFT formalism for excitation energies has become very popular for practical 
applications, and is nowadays implemented as an integral part of several widely used soft- 
ware packages in theoretical chemistry— Overall, TDDFT affords a unique balance between 
accuracy and efficiency, allowing the user to study systems that would be impossible to treat 
with traditional wavefunction methods, for example in all-electron studies of the photochem- 
istry of large biomolecules.— ^ To describe localized and delocalized electronic excitations in 
large conjugated and aggregated molecules, TDDFT has also been used along semiempirical 
Hamiltonian models or approaches based on the Kohn-Sham density matrix such as the 
collective electronic oscillator (CEO) model.— ^ 

The broad spectrum of applications of TDDFT for excited states has been recently re- 
viewed by Elliott et al— From the large body of available literature, the following trends for 
molecules have emerged. Transition frequencies, calculated with standard gradient-corrected 
XC functionals are typically good to within 0.4 eV. Excited-states structural properties such 
as bond lengths, vibrational frequencies and dipole moments are essentially as good as those 
of ground-state calculations (about 1 % for bond lengths, and about 5 % for dipoles and 
vibrational frequencies). Most importantly for large systems, the computational costs scale 
as iV 3 , versus iV 5 for wavefunction methods of comparable accuracy (eg CCSD, CASSCF). 

While there exist efficient iterative algorithms for solving the full eigenvalue problem ([1]), 
it is nevertheless useful to consider approximations, since these can lead to further insight. 
The Tamm-Dancoff approximation (TDA) neglects the off-diagonal matrices K in Equation 
(CQ), which results in the simpler eigenvalue problem 



Physically, the TDA is the TDDFT analog of the configuration interaction singles (CIS) 



a' a' Kiacr,i' a' a' £i 



a'i'a' = 




(4) 



AX — u^tdaX . 



(5) 
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method.— The TDA has some technical advantages away from ground-state equilibrium 
geometry, as discussed by Casida et al. in Ref. Q. 

In an even more drastic approximation, Equation (j3J) is truncated down to a 1 x 1 matrix. 
This yields the so-called small-matrix approximation (SMA),— which for spin-saturated 
systems is given by 

W SMA = <*& + ^iaK ia ,ia- (6) 

This can be further approximated if the correction to the bare Kohn-Sham excitation energy 
u) iarj is small, which is known as single-pole approximation (SPA):— ^ 

i^spa = Wi a + 2K ia)ia . (7) 

This approximation can also be viewed as a TDA for a two-level system. Below, we shall 
present an alternative, more direct derivation of the SMA and SPA. It turns out that these 
two approximations perform surprisingly well for simple closed-shell atomic systems.— This 
might perhaps be viewed as merely a curiosity — after all, one can describe such systems 
with the full Casida TDDFT formalism without resorting to any approximation. However, 
the SMA and SPA can prove very useful for situations such as extended systems where the 
full Casida formalism cannot be easily applied, and we shall hence focus on such cases. 

The purpose of this paper is to present a discussion of simplified TDDFT approaches 
to excitation energies in bulk materials and quantum wells. We will begin with a simple 
and pedagogical derivation of the small-matrix and single-pole approximation for a two- 
level Kohn-Sham system. We will then show how these expressions can be easily modified 
and generalized for the case of collective excitations in extended systems. In this way, one 
can arrive at a simple treatment of excitonic effects in bulk insulators, and plasmon-like 
excitations in doped semiconductor nanostructures. Some examples will be discussed. 

II. EXCITATION ENERGIES OF A TWO-LEVEL KOHN-SHAM SYSTEM 

Let us consider a two- level system consisting of two orbitals y?i(r) and tp<i (r) which are 
eigenstates of the static Kohn-Sham Hamiltonian 

H° = + Kxt(r) + / d 3 r^± + V xc [n ](r) , (8) 

where V xc is the exchange-correlation (xc) potential, a functional of the ground-state density 
Hq(t). We assume that <pi is doubly occupied and y? 2 is empty. 
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Now consider a weak perturbation XH'(r,t) acting on the system. According to time- 
dependent perturbation theory, the time evolution of the system is given by 

<p(r, t) = ci(t)v?i(r) + \c 2 {t)ip 2 {r) . (9) 

Let us construct the density matrix of the system as follows: 



P(t) 



^ Pll Api2 \ I |ci| 2 \c\c 2 



(10) 



V 



Ap 2 i A 2 p 2 2 J \ AciCg A 2 |c 2 | 2 
where we explicitly indicate the order of the perturbation through orders of A. The density 
matrix obeys the following equation of motion: 

p = -i[H,p} = -i[H° + \H'(t),p). (11) 

Dropping terms of order A 2 , this yields the time evolution of the off-diagonal elements of 
the density matrix as follows: 

p 12 = -i[(H° n - H° 22 )p 12 - H[ 2 p n ] (12) 
p 21 = i{{H° n - H° 22 )p 21 - H' 2lPll ], (13) 

where = J d 3 ripi(r)H°(pi(r) and similar for all other matrix elements of H° and H'. 
Since p n = 1 + 0(X 2 ), and defining H 22 — H® x = uj 2 \ (the bare Kohn-Sham excitation 
energy), this simplifies to 

f>12 = A^2\P\2 + H' 12 ] (14) 

P21 = -i[u 21 p 21 + H 21 ]. (15) 

Next, we make the ansatz (which will be justified later) 

p l2 (t) = ~p X2 {uj)e~ wt + ~p 12 {-u)e^ , (16) 

and similar for P21, H[ 2 , and H' 21 . This gives 

-up l2 {u) = [uj 21 p 12 (uj) + H' l2 (u)] (17) 
-up 21 (uj) = -\u 21 p 21 {u) + H' 21 {u)\ (18) 

and an additional two equation for p\ 2 {—uj) and p 2 \(—uj) which do not contain any new 
information. Adding Equations (TT71) and (ITS]) gives 

p 12 {LO) + p 21 {UJ) = ; . (19) 

Co>21 + OJ UJ 2 \ — UJ 



Let us now consider the perturbing Hamiltonian: 

1 



H'(t,u) 



d 3 r' 



r — r 



- + / xc (r, r',u) 



Sn(r', to) , 



(20) 



where 

5n(r,u) = 2(p 1 (r)(p 2 (r)[p 12 {uj) + p2i{u)}- (21) 

Notice that we do not consider an external perturbation, only the linearized Hartree and XC 
potentials. We are thus looking for an "eigenmode" of the system in a steady-state. This 
justifies the ansatz ffTBT) made above. We define the double matrix element 

1 



F(u) = 2 J d 3 r J d 3 r' Vl (r) V2 (r) — ^— + / xt (r. r'. 
and Equation (|T9|) becomes 

Px 2 {ijo) + P2i(u) = -F{u)\p 12 {ijo) + p2i{u)\ 
Canceling pi 2 + p 2 i on both sides results in 



(22) 



1 1 

+ 



(23) 



2uj 



21 



^21 - ^ 2 



F(u) 



which gives the final result 



uj 2 = lo 12 + 2u 2 iF(uj) 



(24) 



(25) 



This is the same as the SMA, Equation ([6]). From the point of view of our two- level system, 
our derivation shows that the SMA considers the excitation 1 —>■ 2 as well as the de-excitation 
2 — > 1. The SPA, Equation (I7j), on the other hand, only includes the excitation 1 — > 2 [it 
is obtained by ignoring the first pole in Equation (1231) ]. In general, the TDA ([5]) ignores all 
de-excitations. 



III. OPTICAL EXCITATIONS IN BULK INSULATORS 

The standard ab-initio treatment of excitation processes in insulators and semiconductors, 
including correlation-induced screening, is based on many-body Green's function techniques 
such as GW and the Bethe-Salpeter equation.— However, TDDFT has recently emerged 
as an alternative, computationally convenient approach to electronic excitation processes 
in materials.— As discussed above, excitation energies can be calculated in principle 
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24. a n approximate 



25 



uses an exact- 



FIG. 1: Schematic illustration of interband transitions in an insulator. The thick arrows represent 
single-particle transitions between the valence (v) and the conduction (c) band. Via dynamical 
many-body effects, these transitions are coupled (thin blue lines). Excitons can thus be viewed as 
collective excitations, where the electron-hole interactions involve the entire bands. 

exactly in TDDFT, provided the XC kernel f xc (r, r',u) is known, and one starts with an 
exact ground-state calculation. This statement is true not only for fin ite systems such 
as molecules, but remains valid for extended systems. In Refs. 
f xc was constructed from many-body Green's functions, whereas Ref. 
exchange approach, including a cutoff in wavevector space which mimics screening of the 
Coulomb interaction.— These studies have established that TDDFT is capable of describing 
excitonic effects in solids, although one has to use XC functionals that go beyond the more 
common ones such as the adiabatic local-density approximation.- The resulting agreement 
with experimental data is excellent,— but the technical effort is not significantly less than 
that of standard many-body approaches. 

In the following, let us for simplicity consider a two-band insulator with valence and 
conduction band energies and e£ and associated Bloch functions Vw( r ) and ^ cq (r). The 
single-particle interband transition energies are simply given by ui^ = e£ — ejj.. A TDDFT 
formalism for nonlinear ultrafast interband excitations has recently been developed,— which 
describes the electron dynamics of a two-band system via a generalization of the semiconduc- 
tor Bloch equations. This work has shown that TDDFT is capable of producing excitonic 
signatures in the interband dynamics. We now briefly outline a simpler and more direct 
TDDFT approach for excitonic binding energies.— 

To describe excitonic effects in the interband absorption spectrum of insulators, one 
needs to go beyond the single-particle picture and include dynamical many-body effects. In 
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TDDFT, this is accomplished through the dynamical XC kernel / xc (r, r', to). Let us illustrate 
this for the simple case of our two-band insulator. The SPA of Equation ([7]) refers to two 
discrete levels, but it can be generalized to two entire bands as follows: 

E [<*«i + F ^)] CqM = "Ck(w) , (26) 
q 

where 

F kq M = ^ J a d 3 r J n d 3 r' ^ k (r)^ k (r)/ xc (r, r', o;)^ q (r')^ q (r') . (27) 

Thus, in contrast to the simple SPA Equation ([7j) which gives an explicit expression for the 
excitation energy, Equation (1261) represents an eigenvalue equation which couples interband 
transitions of different wavevectors via the coupling matrix Fk q , Equation ( 1271) . This shows 
explicitly that excitons are collective excitations which involve not only the states at the 
valence-band maximum and conduction-band minimum (for direct excitons), but the states 
of the entire bands. A schematic illustration is given in Figure HJ Explicit evaluations of 
Equations (|26l) and (1271) are currently in progress.— 



The TDDFT formalism of Ref. 



28l can in principle also be used to treat higher-order 



correlation effects such as biexciton formation. These effects are especially important in 
the description of the ultrafast nonlinear optical response in semiconductors (bulk and 
nanostructured).— Such nonlinear effects, which are entirely governed by dynamical cor- 
relations, can be expected to put severe demands on the XC functional, and will be the 
subject of future studies. 



IV. INTERSUBBAND PLASMONS IN QUANTUM WELLS 

The terahertz frequency regime is scientifically rich, but despite recent progress it is 
technologically still underdeveloped.— 1 Subband spacings in typical III-V quantum wells are 
of the order of a few tens of meV. Since a photon energy of 10 meV corresponds to a frequency 
of 2.4 THz, intersubband transitions in quantum wells appear as natural candidates for 
device applications to fill the "Terahertz gap" in the electromagnetic spectrum. A schematic 
illustration of intersubband transitions is given in the top panel of figure [2j 

The fundamental coupling mechanism between electromagnetic waves and carriers in 
a doped semiconductor quantum well is through a collective excitation, the so-called in- 
tersubband plasmon. To date, most applications of TDDFT linear response theory in 
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FIG. 2: Top: Schematic illustration of interband transitions (blue arrow) and intersubband transi- 
tions (red arrows) in a square quantum well. Bottom: Frequency versus in-plane wave vector k\\ of 
intersubband excitations in a quantum well. Collective charge and spin plasmon excitations (CP 
and SP) are separated from the region of incoherent single-particle excitations. 



quantum wells have been concerned with calculating intersubband plasmon dispersions and 
linewidths.— &^£k£!L£Ii%LM. Due to its highly collective nature, and since the dynamics is 
essentially one-dimensional, intersubband plasmons in quantum wells are an ideal testing 
ground for new TDDFT approaches. A review of recent applications of TDDFT to the 
electron dynamics in semiconductor nanostructure is given in Ref. |42. 

Let us now extend our discussion to include charge as well as spin plasmon excitations. 
The physical picture of these excitations is that in a charge plasmon the spin-up and spin- 
down electrons in the quantum well move collectively in phase, whereas in a spin plasmon 
they move 180 degrees out of phase. These excitations have a dependence on the in-plane 
wavevector k||, which is schematically illustrated in Fig. [2j In the u — k\\ plane, one 
distinguishes an area of incoherent single-particle excitations (also known as Landau damp- 
ing), and two discrete plasmon branches. The charge plasmon lies above the single-particle 
regime, and the spin plasmon lies below. Experimentally, this is confirmed via inelastic light 



scattering spectroscopy. 



43 
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We shall in the following discuss how TDDFT describes the charge and spin plasmons 
in quantum wells. We adopt the effective-mass approximation for electronic states in doped 
semiconductor nanostructures. Let N s be the electronic sheet density (the number of elec- 
trons per area), which is chosen so that only the lowest subband is occupied. It is assumed 
that the direction of growth, and thus of quantum confinement, of the quantum well is along 
the z-axis, and electrons are free to move in the x-y plane. For simplicity, we consider only 
the first and the second subband. 

The intersubband charge and spin plasmon frequencies u CtS are determined by the condi- 
tion that the following 2x2 matrix equation, 



^M ff(J /(k||,^)/v(k||,^) = Xh a (k\\,w) , 



(28) 



has eigenvalue A = 1. The matrix elements are given by 

"2tt 



M aa , (q, u) = F 21 (q, u)ljdzj dz' 



- At / dz 

"21 



4 \ / P(z) 

3 J \ n{z) 



2 



where 



P(z) = ifi{z)V z ip 2 {z) - (p 2 {z)V z (fi{z) , 



and 



F21 (q,^) 



d 2 k 



f{e x + k 2 /2) 



+ 



(29) 



(30) 



(31) 



(2tt) 2 [u + ^-d-k+^ + iS -to + lu 21 - q • k + f - iS J ' 

Here, (pi(z) and <^ 2 (z) are the envelope functions of the first two electronic subbands, e\ t2 
are the associated Kohn-Sham energy eigenvalues, and u 2 \ = e 2 — E\ is the bare subband 
spacing (assuming parabolic subbands). A detailed derivation of Equations (l29]) - (l3Tj) will 
be published elsewhere.— 

The matrix elements M aa i of Equation (1291) contain various XC contributions. The first 
XC term features f^^~(z, z') in the adiabatic local-density approximation (ALDA).— This 
term is frequency independent and real, and therefore by itself gives rise to plasmon excita- 
tions that have no imaginary part, and thus have an infinite lifetime. 
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In reality, since quantum wells are extended systems, plasmons are intrinsically damped. 
Finite plasmon linewidths arise from the other XC terms in Equation (1291) . featuring the 
viscosity coefficients T] XCta(T i and ( XC) aa> and the spin Coulomb drag transresistivity p-^. These 
terms originate^ 1 ^ 1 ^ in the context of time-dependent current-DFT, and have been exten- 
sively studied for quantum wells.— is&aSjilL 

One can simplify things for the case of vanishing in-plane wavevector, ku =0, and ignoring 
the XC contributions beyond ALDA (which in general have only a minor influence on the 
plasmon frequency). In that case, one obtains a straightforward generalization of the SMA 
for quantum wells as follows: 

u 2 = u 2 2l + 2u 21 {S^±S n ) , (32) 



where 



S ™' = T J dZ i dZ 'y ~ ^ " A + f£°*( z > z')\Vi{z) V2 {z) Vl {zi) V2 { Z i) . (33) 
The Hartree contribution in S a(T > is known as "depolarization shift" , and the XC contribution 
is sometimes (somewhat misleadingly) called "excitonic shift". The Hartree part always 
induces an up-shift of the plasmon frequency with respect to u 2 i, and the XC part gives 
a smaller down-shift. In the charge plasmon, the positive shift dominates, but for the 
spin plasmon, the Hartree parts cancel out in Eq. fl32|) and the spin plasmon frequency 
is redshifted with respect to the single-particle excitation region, see Fig. [2j This is a 
remarkable result: the existence of the collective intersubband spin plasmon is purely a 
consequence of XC effects. 

The picture that emerges from these studies is that TDDFT has been extremely success- 
ful in describing excitations in doped semiconductor nanostructures (which are essentially 
metallic systems, in contrast with the insulators discussed in the previous section), already 
at the level of the simplest XC approximations. Currently, work is in progress where the 
formalism is generalized to more realistic descriptions of the semiconductor materials, for 
example including spin-orbit coupling effects.— 



V. CONCLUSION AND CHALLENGES FOR TDDFT 



In this paper, we have given a brief overview of the TDDFT methodology to calculate 
excitation energies in materials. The basic formalism, which is widely used in popular 
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quantum chemistry software packages, is based on Casida's equation (pQ). This approach is 
universal and in principle capable of producing the exact excitation energies of any system 
and material, provided that XC effects are properly accounted for. 

Besides the full formalism, it is instructive to consider simplified approaches which are 
based, in one way or another, on truncations of the matrices A and K of Equation (pQ), 
leading to the SMA or SPA. We have shown how these approaches can be formulated for 
excitations in bulk materials and doped semiconductor quantum wells. This allows for a 
straightforward discussion of excitonic effects in insulators, and intersubband plasmons in 
quantum wells. In both cases we have given explicit expressions for calculating the collective 
excitations. For extended systems, TDDFT not only gives the excitation energies in principle 
exactly, but also the linewidths. 

However, the success of TDDFT critically depends upon the quality of the XC function- 
al used. Simple local and semilocal functionals, such as the LDA or GGA, usually work 
fine for most applications for molecular systems.— On the other hand, there are situations 
which are difficult to deal with, namely a proper description of double excitations or of 
charge-transfer excitations in large molecules. In these cases, one needs to go beyond the 
adiabatic approximation and use an XC kernel f xc (r, r', u) which is explicitly frequency- 
dependent.— 1 ^ The exploration of non-adiabatic XC effects in TDDFT remains a subject of 
intense investigation.— 

For extended insulating or semiconducting systems (inorganic as well as organic), the 
situation is even more complicated. A proper description of excitonic effects in the optical 
absorption of such materials requires long-range XC kernels.— Standard XC functionals 
such as LDA and GGA fail to produce any excitonic effects, and are in fact completely 
incapable of shifting the Kohn-Sham band edge.— >^ It has been demonstrated2^2£>2&2£ that 
one can construct XC kernels which produce absorption spectra in excellent agreement 
with experiment. However, these kernels are technically involved and do not lead to major 
computational advantages compared to Green's function based techniques. It remains an 
important task for TDDFT to construct parameter-free expressions for f xc which accurately 
reproduce excitonic effects in solids, yet are simple enough to implement in practice. Work 
along these lines is in progress. 
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